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A detailed study is presented of the counterrotating model (CRM) for generic electrovacuum static 
axially symmetric relativistic thin disks without radial pressure. We find a general constraint over the 
counterrotating tangential velocities needed to cast the surface energy-momentum tensor of the disk 
as the superposition of two counterrotating charged dust fluids. We also find explicit expressions for 
the energy densities, charge densities and velocities of the counterrotating fluids. We then show that 
this constraint can be satisfied if we take the two counterrotating streams as circulating along electro- 
geodesies. However, we show that, in general, it is not possible to take the two counterrotating fluids 
as circulating along electro-geodesies nor take the two counterrotating tangential velocities as equal 
and opposite. Four simple families of models of counterrotating charged disks based on Chazy- 
Curzon-like, Zipoy-Voorhees-like, Bonnor-Sackfield-like and Kerr-like electrovacuum solutions are 
OO ' considered where we obtain some disks with a CRM well behaved. The models are constructed 

y^ [ using the well-known "displace, cut and reflect" method extended to solutions of vacuum Einstein- 

^>— t ■ Maxwell equations. 
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I. INTRODUCTION 



O < Stationary or static axially symmetric exact solutions of Einstein equations describing relativistic thin disks are 

5^ ' of great astrophysical importance since can be used as models of certain stars, galaxies and accretion disks. These 

{^, were first studied by Bonnor and Sackfield [l|, obtainingpressureless static disks, and then by Morgan and Morgan, 

obtaining static disks with and without radial pressure @, Q- I^i connection with gravitational collapse, disks were 

first studied by Chamorro, Gregory and Stewart [4|. Disks with radial tension have been also studied [9]. Several 

^ classes of exact solutions of the Einstein field equations corresponding to static and stationary thin disks have been 

CN obtained by different authors 0, 0, i, i, [M [H III, [H, 0, d, llll , with or without radial pressure. 

0^ ^ In the case of static disks without radial pressure, there are two common interpretations. The stability of these 

l^ ' models can be explained by either assuming the existence of hoop stresses or that the particles on the disk plane 

^^ , move under the action of their own gravitational field in such a way that as many particles move clockwise as 

C ' counterclockwise. This last interpretation, the "counterrotating model" (CRM), is frequently made since it can 

C^ , be invoked to mimic true rotational effects. Even though this interpretation can be seen as a device, there are 

observational evidence of disks made of streams of rotating and counterrotating matter [ij, 13 ■ 

Disk sources for stationary axially symmetric spacetimes with magnetic fields are also of astrophysical importance 
mainly in the study of neutron stars, white dwarfs and galaxy formation. Although disks with electric fields do not 
have clear astrophysical importance, their study may be of interes in the context of exact solutions. Thin disk have 
been discussed as sources for Kerr- Newman fields |19l] , magnetostatic axisymmetric fields .20.] and conformastationary 
cd ' metrics [21| . Following the Ref. [l3| the resultating disks can be interpreted either as rings with internal pressure and 
currents or as two counterrotating streams of freely moving charged particles, i.e. which move along electro-geodesies 
(solution to the geodesic equation in the presence of a Lorentz force) . 

In all the above cases, the disks are obtained by an "inverse problem" approach, called by Synge the "^-met/iod" [23 |. 
The method works as follows: a solution of the vacuum Einstein equations is taken, such that there is a discontinuity 
in the derivatives of the metric tensor on the plane of the disk, and the energy-momentum tensor is obtained from the 
Einstein equations. The physical properties of the matter distribution are then studied by an analysis of the surface 
energy-momentum tensor so obtained. On the other hand, a "direct problem" approach, called by Synge the "T- 
method^ , is also used by other authors [23, [2I, [23, [2a, [23, [23, [3 by taking a given surface energy-momentum tensor 
and solving the Einstein equations in the matter region. The inner solution is then used to obtain boundary data for 
the vacuum field equations in the outer region. The T-mcthod is used in Re f. .23.1 t o obtain the solution to the problem 
of a (one-component) uniformly rotating disk of dust, and in Refs. [2J, [2a, [2a [23, [23, [231 ^^ generate counterrotating 
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dust disks, but no condition is imposed there about the (electro-) geodesic motion of the two counterrotating streams. 

The aim of this paper is to perform a detailed study of the CRM for generic electrovacuum static axially symmetric 
relativistic thin disks without radial pressure. The counterrotating model for the case of static thin disks without 
electric or magnetic fields was recently studied in [30| , so the material presented here is a continuation of the mentioned 
work. The paper is organized as follows. In Sec. II we present a summary of the procedure to obtain thin disks models 
with a purely azimuthal pressure and currents, using the well-known "displace, cut and reflect" method extended to 
solutions of Einstein-Maxwell equations. In particular, we obtain expressions for the surface energy-momentum tensor 
and the current density of the disks. 

Next, in Sec. Ill, the disks are interpreted in terms of the CRM. We find a general constraint over the counterrotating 
tangential velocities needed to cast the surface energy-momentum tensor of the disk as the superposition of two 
counterrotating charged dust fluids. We also find explicit expressions for the energy densities, current densities 
and velocities of the counterrotating fluids. We then show that this constraint can be satisfied if we take the two 
counterrotating streams as circulating along electro-geodesies. However, we show that, in general, it is not possible to 
take the two counterrotating fluids as circulating along electro-geodesies nor take the two counterrotating tangential 
velocities as equal and opposite. 

In the following section. Sec. IV, four simple families of models of counterrotating charged disks based on Chazy- 
Curzon-like, Zipoy-Voorhees-like, Bonnor-Sackfield-like and Kerr-like metrics are presented where we obtain some 
disks with a CRM well behaved. In particular, we study the tangential velocities, mass and electric charge densities of 
both streams in the pure static and electrostatic (magnetostatic) cases. Also the stability against radial perturbation 
is analyzed. Finally, in Sec. V, we summarize our main results. 

II. ELECTROVACUUM STATIC RELATIVISTIC DISKS 

In this section we present a summary of the procedure to obtain electrovacuum static axially symmetric relativistic 
thin disks. The simplest metric to describe a static axially symmetric spacetime is the Weyl's line element 

ds' ^ - e^'-dt' + e-'-[r'dip' + e^\dr^ + dz% (1) 

where v and A are functions of r and z only. The vacuum Einstein-Maxwell equations, in geometrized units such that 
SttG = c = /io = Eq = I7 are given by 

Rab = Tab, (2a) 

Tab = FacF.'-^gabFcdF^'', (2b) 

F«% = 0, (2c) 

Fab = A,a - Aa,b, (2d) 

where all symbols are understood. 

For the metric ([T]), the Einstein-Maxwell equations in vacuum are equivalent the complex Ernst equations [3l| 

/A£ = (V£ + 2$*V$) • V£, (3a) 

/A$ == (V£ + 2<I>*V$) • V$, (3b) 

with £ = £* (static spacetime), where A and V are the standard differential operators in cylindrical coordinates and 
/ = e^'^. The metric functions are obtained via 

/ = £ + $$*, (4a) 

^.C = ^(^,C + 2<i>*<i>,c)(f.C + 2$$:c)- ^'^X'^h (4b) 



where a/2C — r + iz^ so that ^/2d,Q = d_r — id^z- The electromagnetic potencials are related to $ via 

At = V2Re$, (5a) 

A^,C = \/2iy(Im$),c. (5b) 

In order to obtain a solution of ^ representing a thin disk at z = 0, we assume that the components of the metric 
tensor are continuous across the disk, but their first derivates discontinuous on the plane z — Q, with discontinuity 
functions 



bab — gabA^^o+ ~ 9ab,z\^^^^^ — '^ Qab,. 



z = 0+' 



Thus, the Einstein-Maxwell equations yield an energy-momentum tensor T^ = Q\ 5{z) and a planar current density 
Jq — ja5{z) = —2Aa^z5{z)^ where 5{z) is the usual Dirac function with support on the disk and 

is the distributional energy- momentum tensor. The "true" surface energy- momentum tensor (SEMT) of the disk, 5^, 
can be obtained through the relation 

S'a = JTU^n = e'-^g^, (6) 

where dsn = ^fgTz dz is the "physical measure" of length in the direction normal to the disk, and the surface current 
density as ja = e^~'^ja- For the metric ([T]), the non-zero components of 5*^ and the current density are 

5° = 26"^-^ {A„ - 2i,,z } , (7a) 

Si = 2e^-^A„, (7b) 

and 

}t = -2e^-^A,,, (8a) 

iv - -'2e''-^A^^z, (8b) 

where all the quantities are evaluated at z = 0+. 

With an orthonormal tetrad ea'' = {V'', W'', X^ Y''}, where 

V" = e-" (1,0,0,0), (9a) 

W'' ^ - (0,1,0,0), (9b) 

r 

X^ = e''-^(0, 0,1,0), (9c) 

ya ^ e''-^(o, 0,0,1), (9d) 

we can write the metric and the SEMT in the canonical forms 

gab = -VaVb + WaWb + XaXb + YaYh, (10a) 

Sab = eVaVb + P^WaWb , (10b) 

where 

e = -5° , p^ ^ Si , (11) 

are, respectively, the energy density and the azimuthal pressure of the disk. 



III. THE COUNTERROTATING MODEL 

We now consider, based on references [32| and [331, the possibility that the SEMT S"''' and the current density 
j° can be written as the superposition of two counterrotating fluids that circulate in opposite directions; that is, we 
assume 

S"'' = Sf + SI'', (12a) 

r = j^+j-, (12b) 

where the quantities in the right-hand side are, respectively, the SEMT and the current density of the prograde and 
retrograde counterrotating fluids. 

Let C/J — (f/^, C/±, 0, 0) be the velocity vectors of the two counterrotating fluids. In order to do the decomposition 
(|12ap and (|12bp we project the velocity vectors onto the tetrad ea'', using the relations [3J| 

C/| = e\Ui , Ul = C/4e^^ (13) 

With the tetrad ^ we can write 

Ui = Yl±^^, (M) 

and thus 



V = 77 77 , (15a) 



V- 


- vlv+ui - y^i - vl\j-U^ 


u+-u_ 


^^ 


-vlu^- ^1-V^_U1 



W^ = ^ --^- , (15b) 



where U± = U^/U±. are the tangential velocities of the fluids with respect to the tetrad. 
Using (Uni), we can write the SEMT as 



S 



,, /(U_,U_)(1-U2_)C/-C/| ^ /(U+,U+)(l-U2_)C/-f/^ 



(U+-U_)2 (U+-U_)2 

/(U+,U_)[(1 - Vl){l - V'i)]i{U^Ut + uiuX) 
(U+-U-)2 ' 

where 

/(Ui,U2) = eUiU2+p^ . (16) 

Clearly, in order to cast the SEMT in the form (|12ap . the mixed term must be absent and therefore the counterro- 
tating tangential velocities must be related by 

/(U+,U_) = 0, (17) 

where we assume that |U±| ^ 1. Then, assuming a given choice for the counterrotating velocities in agreement with 
the above relation, we can write the SEMT as (|12ap with 

Sf = e± U^Ui, (18) 

so that we have two counterrotating dust fluids with energy densities given by 



e± = 



1-U^ 



U^-U± 



U^e. (19) 



Thus the SEMT S"'^ can be written as the superposition of two counterrotating dust streams if, and only if, the 
constraint pT]) admits a solution such that U+ ^ U_. This result is completely equivalent to the necessary and 
sufficient condition obtained in reference [33| ■ 

Similarly, we can write the current density as (J12bp with 



J± =<T±U±, 

where a± are the counterrotating electric charge densities of the fluids which are given by 



cr± 



1 



ul 



U+-U, 



J 



J'^ 



lyi yo 



U^). 



(20) 



(21) 



Another quantity related with the counterrotating motion is the specific angular momentum of a particle rotating 
at a radius r, defined as h± = gcpipU^. We can write 



h+ = 



re-''U± 



1 



VI 



(22) 



This quantity can be used to analyze the stability of the disks against radial perturbations. The condition of stability, 



dr 



> 0, 



(23) 



is an extension of the Raylcigh criteria of stability of a fluid in rest in a gravitational field [35| . 

As we can see from Eqs. (fT4l) . ([19]), (|2T|) and ([22]) . all the physical quantities asociated with the CRM depend 
of the counterrotating tangential velocities U±. However, the constraint (fTT]) does not determine U± uniquely, and 
so we need to imposse some additional requeriment in order to obtain a complete determination of the tangential 
velocities, leading so to a well defined CRM. A possibility, commonly assumed, is to take the two counterrotating 
fluids as circulating along electro-geodesies 



-e±9ab,rUlU'^ = -a±FraU^. 



(24) 



The conservation laws (the Bianchi identities) at the disk then imply that the two counterrotating streams interact 
only by means of gravitational and electromagnetic forces. 

Let u>± = Uj./Uj. be the angular velocities of the particles. In terms of uj± we get 



U± 



and so, using dH]), ^ and ^, ^ takes the form 



so that 



V° 



W^ 



^±, 



(25) 



iO^{gil,riO^ +.900,r) 



UJ+ 






+ UJ- 


2 

e 




UJ^LO^ 


= - 





ejl^t.r +P^]oA^,r 



<^gOO,r + 2jo^t,r 



Pv 
e 



(26a) 



(26b) 



where we have used the Einstein-Maxwell equation (|4b|) and the expressions (|7al) - (|8b[) for the SEMT and the current 
density. From the last expression follows immediately that /(U+,U_) vanishes and so the electro-geodesic velocities 
agree with (fT7)l and we have a well defined CRM. 

Note that in general the two fluids circulate with different velocities. However, when the spacetime is electrostatic 
(or magnetostatic) the two electro-geodesic fluids circulate with equal and opposite velocities, so that the constraint 
P?)) is equivalent to 



U^ 



e 



(27) 



as is commonly assum.ed in the works concerning counterrotating disks. We now have two counterrotating charged 
dust streams with equal energy densities 



e± = ^, (28) 



specific angular momenta 



and electric charge densities 



h±=±re-\l^^, (29) 



1 



_ Pv 

1 .. rr 



Oe± =-0^ "JOA/I -■, (30a) 



.™. =±^e^jr^--l, (30b) 

where (Te± and arn± are the electric charge densities in the electrostatic and magnetostatic cases, respectively. The 
velocities are given by (|?7|) . 

Although the electro-geodesic choice leads to a well defined CRM, sometimes the corresponding physical quantities 
may have unphysical behavior. For instance, the expressions (1211) and (j22p may take imaginary values. So we need 
to consider other solutions of Eq. (fT7|) difFerents of the electro-geodesic velocities. Another possibility, commonly 
considered, is to take the two counterrotating fluids not circulating along electro-geodesies but with equal and opposite 
tangential velocities, 

U± = ±U. (31) 

This choice, that imply the existence of additional interactions between the two streams (e.g. colUisions), leads to a 
complete determination of the velocity vectors in such away that the expressions for the velocities, energy densities 
and specific angular momenta of both streams coincide with the electrostatic (magnetostatic) case, but the charge 
densities are differents. However, this can be made only when < |P(^/(t| < 1. If this is not the case, we cannot take 
the two velocities as equal and opposite. In the general case, the two counterrotating streams circulate with different 
velocities and we can write ([T71) as 

U+U_ = -^. (32) 

However, this relation does not determine completely the tangential velocities, and so the CRM is undetermined. 

IV. SOME SIMPLE EXAMPLES OF COUNTERROTATING CHARGED DUST DISKS 

A. CRM for Chazy-Curzon-like disks 



The first family of solutions considered is a Chazy-Curzon-like solution which is given by 

2 



(l + a)e''/P + (1 -a)e-T/p' 



7^ sin^ ' 



(33a) 



A = -^^^, (33b) 

* " (l + a)eV'' + (l-a)e-Wp' ^ ^ 

A^ = V27gcos6i, (33d) 

where a^ = l-f fe^, with \? = p^ + q^, and 7 is a real constant. Here p and q are the electric and magnetic parameters, 
respectively, p and 6 are the spherical coordinates, related to the Weyl coordinates by 

r = psin6', z + zo — pcoa9. (34) 



Note that we have displaced the origin of the z axis in zq. This solution can be generated, in these coordinates, using 
the well-known complex potencial formalism proposed by Ernst 31] from the Chazy-Curzon vacuum solution |36l.l37l|. 
by choosing the parameter q of Ref. [SlJ as complex. For q = we have an electrostatic solution and for p = one 
obtains its magnetostatic analogue [38|. 

From the above expressions we can compute the physical quantities associated with the disks. We obtain 

_ 87Cos6'e"^[(p-7sin2 0)(l + a)eT/''- (p + 7sin^6l)(l-a)e-''/''] 
' " p3[(l + a)e7/p + (l-a)e-7/p]2 ' ^^^^^ 

_ 87^sin^gcos0e-^ 



16V27pcos( 



-A 



•'* (l + a)p2[(l + a)e7/p + (l-a)e-7/p]3' (^^*^) 

• - 8^279 sin^ge-^ 

^■^ " p[(l + a)e7/p + (l-a)e-7/p]- ^''^''^ 

In order to study the behavior of these quantities we perform a graphical analysis of them for disks with 7 = 1, 
zq = 1.5 and p = g = 0, 0.5, 0.1, and 1.5. In Fig. [2a) we show the energy density e (upper curves) and the azimuthal 
pressure 'p^ (lower curves), as functions of r. We see that the energy density presents a maximum at r = and then 
decreases rapidly with r. We also see that the presence of electromagnetic field decreases the energy density at the 
central region of the disk and later increases it. We can observe that the pressure increases rapidly as one moves away 
from the disk center, reaches a maximum and later decreases rapidly. We also observe that the electromagnetic field 
decreases the pressure everywhere on the disk. Next, the charge and electric current densities jt and j,p are represented 
in Fig. [IJ6). jt has a maximum at the disk center and then falls to zero at infinity, whereas j,^ exhibits a similar 
behavior to the pressure. We also cumputed the functions e, "p^^ jt and \^p for other values of the parameters and, in 
all the cases, we found the same behavior. 

We now consider the CRM for the same values of the parameters. All the significant quantities can also be expresed 
in analytic form from the above expressions but the results are so cumbersome that it is best just to analyze them 
graphically. Wc first consider the two counterrotating streams circulating along electro-geodesies. In Fig. [D^a) we plot 
the tangential velocity curves of the counterrotating streams, U+, U_ for disks with 7 = 1, zq = 1.5 and p — q = 0.5, 
0.1, and 1.5, and U^ for 7 = 1, zq = 1.5 and b = 0.5, 0.1, and 1.5. We can see that U+ and U_ increase initially 
and then U+ falls to zero at infinity and always is less than the light velocity, whereas U_ increases monotonously. 
Therefore these disk models are well behaved only at the central regions. However, when the spacetime is electrostatic 
(or magnetostatic) one finds that the velocity U^ for these disks is always less than the light velocity, but the disks 
with b = 0.5 and zq < 0.76 cannot be built from the CRM because U^ > 1 (not shown in the figure). One also finds 
that the inclusion of electric (magnetic) field and the increasing zq make less relativistic these disks. 

In Fig. [2^6) we have drawn the specific angular momenta of counterrotating fluids /i^, h^ for disks with 7 = 1, 
Zq = 1.5 and p = q = 0.5, 0.1, and 1.5, and h^ for 7 = 1, zq = 1-5, b = 0.5, 1.0, 1.5 and 4. We see that there is a 
strong change in the slope of h^ at certain value of r, which means that there is a strong instability there. We also 
find regions with negative slope where the CRM is also unstable. Also, /i?_ presents instability after certain value of 
r. Therefore these disks models are stable only at the central regions. For electrostatic (or magnetostatic) fields the 
specific angular momenta h^ is an increasing monotonous function of r what corresponds to a stable CRM for the 
disks. However, the CRM cannot be applied for b = 4.0 (bottom curve). Thus the inclusion of electric (magnetic) 
field can make unstable these disks against radial perturbations. 

In Fig. [3ja) the plots of the mass densities e^, e_ are shown for disks with 7 = 1, zq = 1.5 and p — q — 0.5, 0.1, 
and 1.5, and e± for 7 = 1, zq = 1.5, b — 0.5, 1.0, and 1.5, as functions of r. The mass density e+ is always positive, 
falling to zero at infinity, whereas e_ becomes negative after some value of r and then falls to zero at infinity. The 
mass densities e± for electrostatic (or magnetostatic) fields are positive everywhere on the disks, falling to zero at 
infinity. We also see that the presence of electromagnetic field decreases the progradee and retrogade mass densities 
at the central region of the disks and later increases them. 

In Fig. [3]J6) the plots of the electric charge densities a+, a- are shown for disks with 7 = 1, zq = 1.5 and p — q — 0.5, 
0.1, and 1.5, and Uei and am± for 7 = 1, zo = 1.5, b — 0.5, 1.0, and 1.5, as functions of r. We see that a+ falls to zero 
at infinity, whereas cr_ becomes later imaginary. ae± and am± present a maximun at the disks center, then decrease 
monotonously, falling to zero at infinity, and are always real quantities. 

As with the electro-geodesic counterrotating fiuids we obtain expressions that becomes imaginary in portions of the 
disks, we also consider non electro-geodesies CRM for these disks. When the two counterrotating streams do not move 
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on electro-geodesies but have equal and opposite tangential velocities the graphics for the velocities, energy densities 
and specific angular momenta of both streams will have the same behavior that in the electrostatic (magnetostatic) 
case. The expressions for the electric charge densities are differents and are plotted in Fig. [H being always real 
quantities. Thus, for Chazy-Curzon-like fields we can build electro-geodesic counterrotating thin disk sources with 
a well behaved central region, whereas when we have counterrotating streams not moving along electro-geodesic 
or electrostatic (magnetostatic) fields we obtain physically acceptable counterrotating disk for many values of the 
parameters. 

B. CRM for Zipoy-Voorhees-like disks 

The second family of solutions considered is a Zipoy-Voorhees-like solution which can be written as 






x^ — 2/2 



(36a) 
(36b) 



{I + a){x + 1)1 + {\-a){x -!)'<' ^^^^' 

A^ = V2kjqy, (36d) 

where a^ = 1 + 6^, with b^ = p^ + q^, and 7 is a real constant. Here p and q are again the electric and magnetic 
parameters, respectively, x and y are the prolate spheroidal coordinates, related to the Weyl coordinates by 

r2 =:fc2(a;2-l)(l-y2)^ z + zo = kxy. (37) 

This solution can also be generated, in these coordinates, using the well-known complex potencial formalism pro- 
posed by Ernst 31] from the Zipoy-Voorhees vacuum solution [39|, |40], also known as the Weyl 7-solution |4ll. |43| . 
by choosing the parameter q of Ref. |3l| as complex. For q = we also have an electrostatic solution and for p = 
one obtains its magnetostatic equivalent [3g|. The case q = and 7 = 1 corresponds to the Reissner-Nordstrom 
solution [44|, in which case a = m/k, p = e/k, with fc^ = m^ — e^, so that a? = I + p'^, m and e being the mass and 
charge parameters, respectively. When p = and 7 = 2 one arrives to a Taub-NUT-like magnetostatic solution. This 
solution can also be generated, in these coordinates, using a well- know theorem proposed by Bonnor (see Ref. |45| ) 
from the Taub-NUT vacuum solution. Note that for 6 = 0, it reduces to the Darmois metric |44| . 

From the above expressions we can compute the physical quantities associated with the disks. We obtain 



e 



8jy{x^ - l)i^'-i\x^ - f)i^''"-^'> 
[(l + a)(.x + l)-r + (l-a)(S-l)')']2 

X {x[{l + a){x + 1)'' + {l-a)(x- iy][x'^ + (7 - l)f - 7] (38a) 

. ^ Sj'xyjl - f){^' - l)h(^'i\x^ ^ f)hii'^^) 

^'^ [(l + a)(S + l)T + (l-a)(x-l)-r] ' ^ ' 

r _ 16727ra(a;' - l)^^(3-7)(^2 _ ^2-)l(7^-2) 

■'* " [(l + a)(.f + l)^ + (l-a)(x-l)7]3 ' ^ ""> 

. _ 4V27g^(l ~ f){x^ - l)^^(i-^)(x2 - y2)|(7--2) 

^^ " [{l + a){x + \)'< + {\-a){x-l)i] ' ^^^""^ 



where e = ke, p^ = kp^ and jt = k}t- x and y are given by 



2x = y/P + {a + l)^ + y/P + {a-iy, (39a) 



2y = y/P + {a + 1)2 - ^r2 + (a-l)2, (39b) 

where f = r/fc and a = zg/fc, with a > 1. 

We consider as an example a particular case, when 7 = 2. We also study Zipoy-Voorhees-like disks for other values 
of 7, but in all cases we found a similar behavior. In Fig. [5] the plots of the quantities e, pep, jt and j,^ are presented 
for disks with a — 2.5 and p = q = 0, 0.5, 1.0, and 1.5, as functions of r. We see that these functions have a similar 
behavior to the previous case. Equally, the relevant quantities of the CRM are shown in following figures for the same 
values of the parameters, also as functions of f. We also first consider the two counterrotating streams circulating 
along electro-geodesies. 

Here the tangential velocity U+ (Fig. EKa)) is also always less than the light velocity, falling to zero at infinity, 
whereas U_ is an increasing monotonously function. Thus these disk models are also just well behaved at the central 
regions. In electrostatic (magnetostatic) case, one also finds that the velocity V"^ for these disks is always less than 
the light velocity, but the disks with b — 0.5 and a < 1.58 cannot be constructed from the CRM because U2 > 1. 
One also finds that the inclusion of electric (magnetic) field and the increasing a make less relativistic these disks. 
Moreover, we see that h'^ and h^ (Fig- HI^)) present instabilities in the same way than the precedent case, and that 
for electrostatic (or magnetostaic) disks (upper curves) the presence of electric (or magnetic) field can make unstable 
the CRM against radial perturbations. Thus the CRM cannot apply for 6 = 4 (bottom curve). 

The mass density e+ (Fig. [7fa)) is also always positive, falling to zero at infinity, whereas e_ becomes negative after 
some value of f and then falls to zero at infinity. We also see that the presence of electromagnetic field decreases the 
progradee and retrogade mass densities at the central region of the disks and later increases them. For electrostatic 
(or magnetostatic) fields e± is positive everywhere on the disks, falling to zero at infinity. We also find that ct+ 
(Fig. ^b)) falls to zero at infinity, whereas (T- becomes later imaginary. The electrostatic and magnetostatic charge 
densities, ae± and a'm±, are always real quantities, falling to zero at infinity. 

Then, when we consider the two counterrotating streams not moving on electro-geodesies but with equal and 
opposite tangential velocities, one also finds that quantities U2, e±, and h^ have the same behavior than in the 
electrostatic (magnetostatic) case. The electric charge densities is shown in Fig. [51 being always real quantities. 
Therefore, for Zipoy-Voorhees-like fields we can build electro-geodesic counterrotating thin disk sources with only a well 
behaved central region, whereas for counterrotating streams not circulating along electro-geodesic or for electrostatic 
(magnetostatic) fields, we obtain disks with a physically acceptable CRM for many values of the parameters. 

C. CRM for Bonnor-Sackfield-like disks 

The third family of solutions considered is a Bonnor-Sackfield-like solution which is given by 

^ (40a) 

(40b) 

(40c) 

(40d) 

where a^ = 1-1-6^, with &2 — p2 ^ ^2^ ^^^ 7 is a real constant. Here p and q are again the electric and magnetic 
parameters, respectively, u and v are the oblate spheroidal coordinates, related to the Weyl coordinates by 

r^ = k^{u^ + l){l-v^), z + Zo = kuv. (41) 

As in the precedent cases this solution can be generated, in these coordinates, following Ernst's method [3l| from the 
Bonnor-Sackfield vacuum solution |l|]. 





(l + a)eT™t-'« + (1 -a)e-T™t"'"' 


A = 


9 


u2 -|-i;2 


5 


At -- 


^p[g7cot-i«_g-7cot-iu] 


(l + a)eT™t"'"-|- (l-a)e-T=°*"'''' 


-ri-ip - 


= \'2kjqv 
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From the above expressions we can compute the physical quantities associated with the disks. We obtain 



87^(^2 + 1) — 



(u + u)t-+2[(1 + a)e'^'=°* '« + (1 -a)e"T^°* '"]2 



Vv 



Jt 



2 

(u + u)T-+2[(l + a)e'>"=°*"'" + (1 -a)e-T™*"'"] 
16-\/27Pw(m + 1)^ 



. _ 4V27gM(l-w2)(u+i)7V2 

(m + u)t-+1[(1 + a)eT™*"' " + (1 - a)e-T™'"'"] 
where e = fee, p^ — kp^p and jt = fcjt. u and w are given by 



V2u = ^[(f2 + a2 _ 1)2 + 4a2]l/2 + (f2 + ^2 - 1), 
\/2u = ^^[(^2 + a2 _ 1)2 + 4q,2]1/2 _ (^2 + ^2 „ 1)^ 



(42a) 

(42b) 

(42c) 
(42d) 

(43a) 
(43b) 



where f — r /k and a = z^jk. 

In Figs. [9lfT2lthe plots of the physical quantities related with the disks are presented for disks with 7 = 1, a = 1.5 
and p = g = 0, 0.5, 1.0, and 1.5, also as functions of f. We see that these functions have the same behavior that the 
previous cases. Therefore, when we consider coimterrotating streams moving along electro-geodesies only the central 
region of the disks is well behaved, whereas in the case of counterrotating streams not circulating along electro- 
geodesies or for electrostatic (magnetostatic) fields we obtain disks with a physically acceptable CRM for many values 
of the parameters. 



D. CRM for Kerr-like disks 

Finally, other family of solutions to the Einstein-Maxwell equations, is a Kerr-like solution which is given by 



At 



A.. 



In 


a2a;2 — 52^2 _ ^ 






(a:r + 1)2 - 62y2 


' 


■^In 


r ? ? 7 2? 1 
ax — by — 1 






[ a2(^2„j^2) 


1 


2V2ra 




(aa;-M)2-62y2' 


V 


'2kq{\ - y^)(ax ^ \) 



a(a2x2 — 62y2 _ 1) 



(44a) 
(44b) 
(44c) 
(44d) 



where c? = 1 -t- ^2 ^ with fe2 = j)2 ^ ^2 ^ Here p and q are also the electric and magnetic parameters, respectively, x and y 
are again the prolate spheroidal coordinates, given by (j37p . For & = this solution also reduces to the Darmois metric. 
Moreover, for g = we have a Kerr-like electrostatic solution. This solution can be generated, in these coordinates, 
using a well- know theorem proposed by Bonnor (see Ref. [43) from the Kerr vacuum solution. And for p = one 
obtains its magnetostatic equivalent. This solution is asymptotically flat and was firts studied by Bonnor '4(?| in a 
different system of coordinates and describes the field of a massive magnetic dipole. 
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The physical quantities associated with the disks can now be written as 



ia^y 



{a^x^ - b'^y'^ - 1)2 [(ax + 1)^ - 62y2]2 

X {(x^ - f)[a{x^ - l)[{ax + if + b^f] ~ 2b'^x{ax + 1)(1 - f)] (45a) 



- 2x{x'' - l)il - f)[{ax + if - b^f]} 
16a'^xy{x^ - 1){1 - y^) 

(a2i2 _ b2y2 _ 1)2 [(a,^ + 1)2 _ 52y2] ' 



(45b) 



r ^ 4V2a4p(x2 - f){x{l - f)[{ax + 1)^ + b^y^] - 2a^2(a^ + i)(^2 „ ^)| 

■"* (a2a:2 - 62^2 _ l)[(a5; + 1)2 _ 52j^2]3 ' ^ ^^) 

. ^ 2V2a^qy{x^ - 1)(1 - y2)(^2 _ ^2)[(^^ ^ i)(3^^^ ^ i) ^ ^2^2] 

•''^ (a2a:2-fe2y2_i)3[(a^ + l)2_52y2] ' ^ ^^^ 

where x and y are given by Eqs. p9ap and (|39b[) . 

In figs. [TSlfTBl the plots of the physical quantities describing the disks are shown for a = 1.7 and p = q = 0, 0.4, 
0.7, and 1.5, as functions of f . The energy density behaves in the opposite way to the previous cases. That is, near 
to the disk center it increases when the electromagnetic field is applied and decreases later. The other quantities 
have a similar behavior to the earlier cases. However, the electric charge density jo after some f takes negative 
values. Likewise, the quantities corresponding to the CRM are shown in following figures for the same values of 
the parameters, also as functions of f. Again, we first consider the two counterrotating streams circulating along 
electro-geodesies. 

Now the tangential velocities U+ and U_ are always less than the light velocity, falling to zero at infinity (Fig. 
I14r a)). Note moreover that the two counterrotating streams initially circulate with different velocities and then 
become equal and opossite. Thus these disks models are well behaved everywhere on the disks. We also found that 
for some values of the parameters cannot be constructed disks from the CRM because U2 > 1. We also note that the 
presence of electromagnetic field can stabilize the CRM against radial perturbations (Fig. [IDJ&)). On the other hand, 
the mass densities e+, e_ and e± (Fig. I15f a)) are also always positive and fall to zero at infinity. 

In Fig. \Wlb) we have plotted the charge densities a+, a-, ae± and am,±- We see that a+ and am± are always 
positive, falling to zero at infinity, whereas a_ can take initially negative values and ae± become negative after some 
value of f and then fall to zero at infinity. Thus, the Kerr-like disk with a = 1.7 and p = q = 1.5 can be constructed 
from the counterrotating streams circulating along electro-geodesies. Finally, the electric charge densities when the 
two streams do not move on electro-geodesies are shown in Fig. 1161 being always real quantities. Therefore, in all the 
cases considered we can build Kerr-like disks with a phisically acceptable CRM for many values of the parameters. 

V. DISCUSSION 

A detailed study of the counterrotating model for generic electrovacuum static axially symmetric relativistic thin 
disks without radial pressure was presented. A general constraint over the counterrotating tangential velocities was 
found, needed to cast the surface energy-momentum tensor of the disk in such a way that it can be interpreted as 
the superposition of two counterrotating charged dust fluids. The constraint found is completely equivalent to the 
necessary and sufficient condition obtained in reference [33]. We next showed that this constraint is satisfied if we 
take the two counterrotating fluids as circulating along electro-geodesies. We note that en general the two fluids 
circulate with different velocities. However, when the spacetime is electrostatic (or magnetostatic) the two geodesic 
fluids circulate with equal and opposite velocities. We also have obtained explicit expressions for the energy densities, 
current densities and velocities of the counterrotating streams in terms of the energy density, azimuthal pressure and 
planar current density of the disks, that are also equivalent to the correspondig expressions in reference 33]. 

Four families of models of counterrotating charged disks were considered in the present work based in simple solutions 
to the vacumm Einsteins-Maxwell equations in the static axisymmetric case generated by conventional solution- 
generating tecniques [44|. When the two counterrotating streams move on electro-geodesisc, Chazy-Curzon-like, 
Zipoy-Voorhees and Bonnor-Sackfleld-like counterrotating disks can be constructed only with a well behaved central 
region, whereas for Kerr-like fields we found for many values of the parameters disks with a physically acceptable 
CRM. On the other hand, when we consider the two counterrotating fluids not circulating on electro-geodesic but with 
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equal and opposite velocities or for electrostatic (or magnetostatic) fields, we found for many values of the parameter 
disks with a well behaved CRM everywhere. In the case of Chazy-Curzon-like, Zipoy-Voorhees and Bonnor-Sackfield- 
like fields, we also saw that the presence of electric (magnetic) field can make unstable against radial perturbations 
the CRM, and conversely, stabilizes the CRM in the case of Kerr-like disks. 

Finally, the generalization of the countcrrotating model presented here to the case with radial pressure is in consid- 
eration. Also, the generalization to rotating thin disks with or without radial pressure in presence of electromagnetic 
fields is being considered. 
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FIG. 1: (a) e (upper curves) and p^ (scaled by a factor of 2.5), as functions of r for Chazy-Curzon-like disks with 7 = 1, 
2o = 1.5 and p = q — (top curves), 0.5, 1.0, and 1.5 (bottom curves), (b) jt (lower curves scaled by a factor of 20) and j^ for 
the same values of the parameters. 
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FIG. 2: (a) U+ (lower curves), U- (upper curves) for Chazy-Curzon-liko disks with 7 = 1, 



zo 



1.5 and p — q 



and 1.5, and U^ for 7 = 1, ^o = 1 and b = 0.5, 1.0, and 1.5. (fe) h\ (scaled by a factor of 1000), h^ (sharp curves 



same values of the parameters and h^ (upper curves scaled by a factor of 1000) for 7 
(bottom curve). 



1, zo = 1.5 and b = 0.5, 
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e+, e-, <;± 
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r 

(a) 



r 

(b) 



FIG. 3: (a) e_ (lower curves), e+ (moved upwards a factor of 0.1) for Chazy-Curzon-likc disks with 7 = 1, zo — 1.5 and 
p — q = 0.5 , 1.0, and 1.5 and e± (upper curves moved upwards a factor of 0.2) for "f = 1, zq — 1.5 and 6 — 0.5, 1.0 and 1.5. 
(6) (7_, 0-+ (curves with crosses), cre± (moved upwards a factor of 0.2) and am± (upper curves moved also upwards a factor of 
0.2) for the same values of the parameters. 




FIG. 4: (7+ (solid curves) and ct_ for not electro-geodesic Chazy-Curzon-like disks with 7 = 1, 20 = 1.5 and p = q — 0.5 (lower 
curves), 1.0, and 1.5 (upper curves). 
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FIG. 5: (a) e (upper curves ) and j5^ (scaled by a factor of 2), as functions of f for Zipoy-Voorhees-like disks with 7 = 2, a = 2.5 
and p — q = (top curves), 0.5, 1.0, and 1.5 (bottom curves). (6) jt (lower curves scaled by a factor of 10) and j^ for the same 
values of the parameters. 
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FIG. 6: (a) U+ (lower curves), U- (upper curves) as functions of f for Zipoy-Voorhees-like disks with j = 2, a — 2.5 and 
p = q = 0.5, 1.0, and 1.5 and U for 7 = 2, a = 2.5 and b — 0.5, 1.0, and 1.5 . (6) h^ (scaled by a factor of 1000), /i_ (sharp 
curves) for the same values of the parameters and n^ (upper curves scaled by a factor of 1000) for 7 = 2, a = 2.5 and b — 0.5, 
1.0, 1.5 and 4.0 (bottom curve). 
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FIG. 7: (a) e_ (lower curves), e+ (moved upwards a factor of 0.1) for Zipoy-Voorhees-like disks with 7 = 2, q = 2.5 and 
p — q = 0.5 , 1.0, and 1.5 and e± (upper curves moved upwards a factor of 0.2) for 7 = 2, a = 2.5 and b — 0.5, 1.0 and 1.5. (6) 
a- (lower curves), a+, d'e± (moved upwards a factor of 0.1) and o"m± (upper curves moved also upwards a factor of 0.1 ) for 
the same values of the parameters. 
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FIG. 8: cT-i- (upper curves) and cr_ for not electro-geodesic Zipoy-Voorhees-like disks with 7 = 2, q = 2.5 and p = q — 0.5 
(lower curves), 1.0, and 1.5 (upper curves). 
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FIG. 9: (a) e (upper curves ) and p,^ (scaled by a factor of 2.5), as functions of f for Bonnor-Sackfield-like disks with 7 = 1, 
a = 1.5 and p = q = (top curves), 0.5, 1.0, and 1.5 (bottom curves), (b) jt (lower curves scaled by a factor of 5) and }^ for 
the same values of the parameters. 
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FIG. 10: (a) U+ (lower curves), U- (upper curves) for Bonnor-Sackficld-like disks with 7 = 1, a = 1.5 and p = q — 0.5, 1.0, 
and 1.5 and U for j — 1, a = 1.5 and b — 0.5, 1.0, and 1.5, as functions of f. (b) h+ (scaled by a factor of 500), /i_ (sharp 
curves) for the same values of the parameters and h? (upper curves scaled by a factor of 400) for ') — 1, a — 1.5 and b = 0.5, 
1.0, 1.5 and 4.0 (bottom curve). 
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FIG. 11; (a) e_ (lower curves), e+ (moved upwards a factor of 0.1) for Bonnor-Sackfield-like disks with 7 = 1, a = 1.5 and 
p — q = 0.5 , 1.0, and 1.5 and e± (upper curves moved upwards a factor of 0.2) for 7 = 1, a = 1.5 and b — 0.5, 1.0 and 1.5. 
(b) cr_ (lower curves), a+, ae± (moved upwards a factor of 0.2) and am± (upper curves moved also upwards a factor of 0.2) for 
the same values of the parameters. 
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FIG. 12: cr+ (upper curves) and cr_ for not electro-geodesic Bonnor-Sackfield-like disks with 7 = 1, a = 1.5 and p — q = 0.5 
(lower curves), 1.0, and 1.5 (upper curves). 
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FIG. 13: (a) e (lower curves ) and p^, as functions of f for Kerr-like disks with a = 1.7 and p = q = (bottom and top curves) 
0.4, 0.7, and 1.5 (top and bottom curves), (fe) jt (upper curves ) and j^ for the same values of the parameters. 
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FIG. 14: (a) U+ and U- (curves with crosses) for Kerr-like disks with a = 1.7 and p = q = 0.4, 0.7, and 1.5, and U'^ for 
a = 1.7 and b — 0.4, 0.7, 1.5, as functions of f . (6) h^ (curves with crosses), h^ for the same values of the parameters and h^ 
(upper curves scaled by a factor of 2) for q = 1.7 and b — 0.4, 0.7, 1.5 and 4.0 (bottom curve). 
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FIG. 15: (a) e_ (lower curves), e+ (moved upwards a factor of 0.1) for Kerr-like disks with a — 1.7 and p — q = 0.4, 0.7, and 
1.5 and e± (upper curves moved upwards a factor of 0.4) for a = 1.7 and b = 0.4, 0.7 and 1.5. (b) cr_ (lower curves), a+, ae± 
(moved upwards a factor of 0.5) and am± (curves with crosses moved also upwards a factor of 0.5) for the same values of the 
parameters. 




FIG. 16: cT-i- (upper curves) and o"- for not electro-geodesic Kerr-like disks with a — 1.7 and p = q — 0.4 (lower curves), 0.7, 
and 1.5 (upper curves). 



